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ABSTRACT 

A model for advection-dominated accretion flows (ADAFs) in a global 
magnetic field is proposed. In contrast to the well known ADAF models in 
which the viscosity of a fluid determines both angular momentum transfer and 
energy dissipation in the flow, the magnetic field and the electric resistivity, 
respectively, control them in this model. A manageable set of analytic solutions 
for the flow and the magnetic field is obtained to vertically non-integrated basic 
equations. This set describes mathematically a fully advective accretion flow 
and, in physically plausible situations for most AGNs, it is also confirmed that 
the radiation cooling estimated on this solution is really negligible compared 
with the internal energy of the flow. 

Subject headings: accretion, accretion disks -- magnetohydrodynamics: MHD 
— galaxies: active — black hole physics 
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1. INTRODUCTION 

In recent years, advection-dominated accretion flows (ADAFs) are drawing much 
attention in astrophysics because the radiation spectrum of Sgr A* is well reproduced in 
the framework of such models (Narayan, Yi & Mahadevan 1995; Manmoto, Mineshige & 
Kusunose 1997; Narayan et al. 1998). The attempts are further made to apply them also 
to some other AGNs and soft X-ray transients (for a review, see Narayan, Mahadevan & 
Quataert 1999). Such a flow may be realized when the rate of mass accretion is smaller 
than the Eddington accretion rate, and the flow becomes optically thin while geometrically 
thick. Almost whole energy dissipated in the flow is advected down to the central object 
owing to an inefficient radiation cooling. 

The basic ideas of such ADAF models has been developed by a number of researchers 
(e.g., Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994, 1995; Abramowicz et al. 1995). 
We call them the "viscous" ADAF models because both angular-momentum transfer and 
energy dissipation in the flow is undertaken by the kinematic viscosity whose size is usually 
specified by the a parameter: v = (2/3)aCsH, where Cs is the sound velocity and H is the 
half-thickness of the flow. Another parameter /3 enters into this model through the assumed 
presence of turbulent magnetic fields (Narayan & Yi 1995): p m = 3(1 — j3) C s 2 p, where p m 
is the turbulent magnetic pressure and p is the mass density. The spectrum of Sgr A* in a 
range from radio up to hard X-ray has been explained by the synchrotron emission from 
high temperature electrons, its inverse- Compton scattered component and bremsstrahlung. 
For the presence of synchrotron emission, the inclusion of turbulent magnetic field is crucial 
in this model. 

Although the above model completely ignores the presence of an ordered magnetic 
field, it is natural to expect the presence of such fields in most AGNs and especially in the 
central region of our Galaxy (e.g., Yusef-Zadeh, Morris & Chance 1984). Taking this fact 
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into account, we propose another ADAF model which shall be called the "resistive" ADAF 
model in order to distinguish it from the viscous ADAF. This name reflects the facts that 
angular momentum is transported by a global magnetic field and energy is released as a 
resistive dissipation of the electric current driven by a rotational motion of the accreting 
plasma. Originally, this model appeared (Kaburaki 1986, 1987) as a magnetic counterpart 
of the standard model (usually called the a-disk model) for optically thick, geometrically 
thin accretion disks (Shakura & Sunyaev 1973). 

Analogously to the a-disk model, this model contains a parameter A which specifies 
the efficiency of dissipation: a~ l oc A 2 , where a~ l is the electric resistivity. The parameter 
A in this model also has a geometrical meaning of the half opening-angle of a disk, which is 
assumed for simplicity to be constant throughout the disk. The set of analytic solutions for 
the flow and the magnetic field was obtained there under the assumption of geometrically 
thin disk (A <C 1), and it describes an accretion flow which has a reduced Keplerian 
rotation and a magnetically confined vertical structure. 

Although the heating process in the disk has been discussed in the above papers, 
the cooling process has not been discussed and implicitly assumed to be in local balance 
with the heating. Since the energy equation is not included explicitly, the solution may be 
inconsistent from a viewpoint of energy transport. In principle, the fraction of advected 
part of the released energy can be calculated from that set of solutions. Actually, however, 
the inclusion of a vertical motion in the model of Kaburaki (1987, hereafter referred to as 
K 87) has made this task rather difficult and obscure. Therefore, in the present paper, we 
intend to obtain a more manageable set of solutions to a similar problem by adding further 
simplifying assumptions of no vertical flows and electric currents. Although such a solution 
has already been obtained in a previous paper (Kaburaki 1986) , it is unsatisfactory because 
the effect of pressure gradient has been omitted without justification. 
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Our new solution obtained in the present paper turns out a posteriori to describe a 
fully advective accretion flow (as a preliminary report, see Kaburaki 1999). In other words, 
the requirements for dynamical balances and the geometry of a flow specify in our model 
also the way of energy transport. The consistency of this solution as an ADAF model in 
most physical situations expected in typical AGNs is also confirmed. Namely, it is shown 
by using this solution that the flow is actually optically thin and the radiation cooling is 
negligible compared with the heating. Therefore, our resistive ADAF model may be useful 
in calculating spectra form various AGNs and in considering stabilities of realistic ADAFs. 
Some results in such applications will appear elsewhere. 

This paper is organized as follows. In §2, a set of resistive MHD equations is first 
introduced as our starting point, and then various assumptions to simplify these equations 
are stated. The resulting equations in spherical polar coordinates are cited in §3, and 
the effects of the individual assumptions and the policy for solving these equations are 
discussed. A set of approximate solutions which we propose as a model of resistive ADAFs 
is written out in §4 and its properties, including energy budget, are examined. §5 is devoted 
to determining the disk edges which are defined conceptually as the limiting radii for the 
validity of this solution. In §6, various quantities predicted by this model are suitably 
scaled to typical situations in most AGNs and the physical consistency of this solution as 
an ADAF model is also checked. Finally, main issues which may be raised about our model 
are discussed in §7. Some remarks concerning a relation between the rotation law and 
geometry of a disk are also described in Appendix. 
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2. BASIC EQUATIONS AND SIMPLIFYING ASSUMPTIONS 

The basic equations adopted in the resistive ADAF model are those of resistive MHD 
which are written in usual notation as follows: 

^ + V-(pv) = 0, (1) 

<9v 11 

_ + (v ■ V)v = — Vp - — [B x (V x B)] + g, (u = 0) (2) 

V-B = 0, (3) 

<9B r 2 

— - = V x (v x B) + AB. (a = const.) (4) 

at Ana 

In the equation of motion (H), the viscosity of the fluid is completely neglected in order to 
make the contrast between the resistive and viscous ADAF models clear. Instead, the term 
of magnetic diffusion is retained in the induction equation (Q) assuming the constancy of 
electrical conductivity a for simplicity. The above set of equations forms an apparently 
closed set in the sense that the number of equations and unknowns are the same. In fact, 
however, the number of equations is insufficient by one. This is because Maxwell's equations 
guarantee that div B = always holds once it is satisfied initially. This point will become 
more clear in §3. Other related quantities of our interest are calculated from the following 
subsidiary equations: the current density, electric field and charge density, from 

1 "1 1 

jVxB, E=i--vxB, q = — V ■ E, (5) 

Air a c An 

respectively. The temperature of the fluid is simply assumed to be common to electrons 
and ions, and is calculated from the ideal gas law 

p= — pT } (6) 

neglecting the radiation pressure (p, is the mean molecular weight and R is the gas constant). 

Hereafter, we employ spherical polar coordinates (r, 8, if) since the gravity is spherically 
symmetric. The main assumptions used in simplifying the above basic equations are, 1) 
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stationarity, 2) axisymmetry, 3) geometrically thin disk, 4) well-developed magnetic disk, 
5) reduced Keplerian rotation and 6) vanishing of vg and jg. All these assumptions except 
the last one are the same as employed in K 87. Although the solution has been obtained 
there without 6), we add it here in order to get a more manageable set of solutions. 

The first two of the above assumptions imply that d/dt = and 8 /dip = in every 
component equation. The most essential assumption is the third one which demands that 
A < 1. This fact allows us to treat A as a smallness parameter in the following discussion. 
In such situations, a derivative with respect to 9 results in a large quantity (introduce 
£ = (9 — 7r/2)/A then 8/ 89 = A _1 <9/<9£). It is worth emphasizing, however, that the thin 
disk condition is satisfied even by rather geometrically thick disks since 1 radian is about 
60° of arc. 

The fourth assumption states that, when the total magnetic field is divided into 
the externally-given seed field B and the component produced in the disk b, we expect 
|b| ^> | Bq | inside the disk except near a disk edge where |b| ~ |B | - The actual 
implication of the fifth assumption is that of large magnetic Reynolds numbers 9ft (strictly 
speaking, 3ft 2 (r) ^> 1) in the disk except near its inner edge. Namely, as it will tern out 
a posteriori from the resulting solutions, the assumption of reduced Keplerian rotation 



{vtp = const, x v k < vk, where vk = JGM/r represents the Kepler velocity) is justified 
when 3ft 2 (r) ^> 1 (see §4). The sixth assumption results in an r-independent accretion rate 
M and the dependence b$ oc r _1 . 

3. LEADING ORDER EQUATIONS IN A 

The set of basic equations is simplified first according to the above assumptions except 
5) and 6). Since the flow is assumed to be geometrically thin, only the leading order terms 



in A are retained in each component of the basic equations. Regarding b r , b v , v r , v v , p and 
p as quantities of order unity in A, we obtain the following component equations, within 
the approximation sin# ~ 1. 

mass continuity: 

d . „ r d 

0. 



d . 2 . r d , 
^ (r pVr) + A dt {fWe 



magnetic flux conservation: 



d 



r dbf> 



Tr(r%) + -r ~^r = 



dr 



A d£ 



equation of motion: 



d vp d 

dr Ar <9£ 



v r — 



1 dp GM 1 
p dr r 2 



A-npr 



^_, d_ ( , 
A 0f v «9r l * 



S + s!fc , + tf 



0, 



dr Ar <9£ 



^ + 



V^Vr 



Airpr 



b r fr^K) 



A di 



induction equation: 



d_ 



v r b e 



47TO-A r <9£ 



9r 



r w r 6 e - 



+ 



1 06/ 



ArraA r <9£ 



9 r / 7 i m ! ^ 

— [r(Vr-UrM] + £ &7 



iV&e 



i^ + 



0, 
1 <%„ 



0. 



(7) 



(9) 
(10) 
(11) 

(12) 
(13) 
(14) 



Ana A r <9£ 

From the equations of mass continuity and flux conservation, it turns out that v$ and bg are 
quantities of order A. Further, the assumption of vq — results in an r-independent radial 
mass flux, r 2 pv r . 



At this stage of approximation, there remain so many terms in the r-component of 
the equation of motion. Only the last term on the the right-hand side can be dropped by 
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a further assumption of jg = (for the component expressions of the current density to 
the leading order in A, see K 87). As confirmed retrospectively, however, p, p oc 3ft 2 while 
v r oc 9ft -1 . Therefore, when 3ft 2 ;^> 1 (the assumption 5), equation (|) reduces to 



GM v l 1 dp 



(15) 



r 2 r p dr 

Owing to the presence of the outward pressure-gradient, the effect of gravity is somewhat 
reduced. Anticipating the same r- and ^-dependences for the gravity and pressure terms, 
we have a reduced Keplerian rotation from this equation. The validity of this anticipation 
is again justified a posteriori from the resulting solution. 

The ^-component of the equation of motion can be integrated to give the "vertical" 

pressure balance for the disk: 

bl 

P + ^ L =P(r), (16) 

where p(r) represent a function of r only. In the above procedure, the magnetic pressure 
due to b r has been omitted in expectation of the relation b v /b r oc 3ft, which is also confirmed 
retrospectively. This equation implies the confinement of the accreting flow into a disk 
structure by the induced toroidal magnetic field which changes sign at the equatorial plane. 
The (^-component of the equation of motion, on the other hand, reduces to 



d(rv v ) 1 b db v 



(17) 



dr 47rA p d£ 

when simplified by the assumption 6). This equation describes the transfer of angular 
momentum by the magnetic stress. The extracted angular momentum is carried away from 
the disk along the poloidal lines of force. 

Both r- and ^-components of the induction equation reduce to 

c 2 1 db r . . 

Vr = ~4^a rb, ae ' (18) 

since the quantities in the brackets of these equations are proportional to the (^-component 
of the electric field E v which should vanish in axisymmetric situations (see Ohm's law in 
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K 87). Thus, the degeneracy of the two equations have become evident. However, the 
degeneracy itself is not a consequence of symmetry but of the nature inherent in Maxwell's 
equations as stated before. We can expect from this equation that v r oc 3ft _1 . 

As for the (^-component of the induction equation, it is also simplified by the 
assumption of Vg — and by the fact that v v is independent of £ since it is proportional to 
the Kepler velocity. Further eliminating be with the aid of magnetic flux conservation, we 
obtain 

In solving this equation, we assume intuitively the proportionality of the first and second 
terms. Afterwards, from the resulting solution, the consistency of such an assumption is 
confirmed. This is the same technique as used in solving equation (|15|). 

4. AN ADAF SOLUTION 

The above set of partial differential equations is incomplete in the sense that the 
number of equations is insufficient by one. Usually, this shortage is supplemented by 
the energy equation or the polytropic one, and it makes the process of solving the set 
of equations rather difficult. Fortunately, however, we can get a solution here without 
introducing another such equation. The set is solved by the method of approximate 
variable separation. Starting from the assumption of b^ oc tanh£, we can obtain almost 
automatically the angular dependences of other quantities according to the equations. The 
only approximation needed is to accept the relation (d 3 /d£ 3 ) tanh£ ~ —2 sech £ which hold 
fairly good around the disk's midplane. Including the method of obtaining the radial parts 
of the unknowns, the process is essentially the same as described in K 87. Therefore, we 
cite below only the final results. 
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The boundary value of the external magnetic field B is taken into account at a 
reference radius r , which may be the inner or outer edge of the accretion disk (they will be 
defined explicitly in the next section). Then the set of solutions for various quantities are 

written as 

/ r \ ~ 3 / 2 
b r (r, = b r (r) sech 2 £ tanh£, bJr) =B [ — ) , (20) 

Vr / 

b e (r, = b e {r) sech 2 £, b e (r) = -B„ {-) , (21) 

b v (r, = -b v (r) tanh£, b v (r) = &(r ) B (-) , (22) 

v r (r, = -firCr) sech 2 e, t5 r (r) = ^ o) (-)"' , (23) 

V3 5R(r ) Vr / 

u*(r, = 0, (24) 

^k(^o) (t \- x I 2 



v<p(r, = iV(r), u v (r) = — -=- ( — J , (25) 

p(r, = p(r) sech 2 £, p(r) = l Q ° ° - , (26) 



Six \r , 
^ v K( r o) W 



p(r, = P(r) sech 2 £, /5(r) = Q _ \*_? ( - ) , (27) 



>(r, = -Jr(r) sech 2 e, j r (r) = -^ ^^ (-) "' , (28) 

47rA r Vr / 

J'e(r, = 0, (29) 

J>, = ~J» sech 4 e, J» = ^ ^ f^)^ , (30) 

^^rMr^iQ' 1 - (31) 

The definition and its r-dependence of the magnetic Reynolds number in our model are 

K Vp ^ s ( r \V2 ^ ^ 7ro-A 2 r wk(^o) 



»(r) = f = ^ = »(r )[- , K(r ) = 7=-f^ ( 32 ) 

6 r v r \r J V3 c 2 

From the latter expression, it turns out that our magnetic Reynolds number is reflecting 
the vertical structure of the disk since Ar and Av k represent typical sizes of the height and 
the velocity in the ^-direction, respectively. 
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Here, we summarize the characteristic features of the above set of solutions. The 
plasma density and pressure are large quantities of order 9? 2 , reflecting the development of 
a disk-like structure of the flow. In order to vertically support this configuration, a large 
toroidal magnetic field b^ of order 9ft is required to appear, and this fact also guarantees 
sufficient extraction of angular momentum from the accreting matter. Controlled by 
this extraction, the plasma falls inwardly. This inflow is understood also as a result of 
resistive diffusion across the poloidal magnetic field lines. Reflecting this fact, v r is inversely 
proportional to 3ft. The presence of a factor 3ft in j r represents that the rotation dominated 
accretion flow acts mainly as a poloidal current generator, from an electrodynamical point 
of view. 

Since the poloidal and toroidal components of vector quantities have different 
r-dependences, our solution is not such a similarity solution as obtained in the framework of 
viscous ADAF (Narayan & Yi 1994). Reflecting the disk's geometry, the rotation velocity is 
a reduced Keplerian (see Appendix for further details) at the points on a surface of constant 
r. However, the plane of rotation is always parallel to the disk's midplane so that the center 
of rotation generally does not coincide with the center of gravitational attraction. It is 
also interesting to see that a polytrope-like relation holds only for the radial functions, i.e., 
p oc p 7 with 7 = 2. The assumption of Vg — assures an r-independent mass accretion rate 

^-^A f^di = y . (33) 

Jo V3 v K (r ) 

A similar thing is assured by the assumption of jg = for the total poloidal current flowing 
in the disk. In order to fulfill the requirement of poloidal current closure, however, non-zero 
jg should be included. Our model, therefore, is shifting this problem to the regions both 
beyond the outer edge and within the inner edge of the disk. 

Based on the above cited solution, we can discuss the energy budget in the accretion 
flow. Since the electric current vector is dominated by its r-component, we have for the 
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local Joule dissipation rate 

o" a 2V3 87r \r / 

where Q = v v /r denotes the angular velocity (within the approximation of sin# ~ 1) and 

a has been eliminated by equation (p2[) . Although the origin of the resistivity need not be 

specified in our model, it is very likely to be of anomalous type since a large current density 

in the disk may cause a plasma turbulence. On the other hand, the advection cooling, 

which is defined as 

I O OX) 

9alv = pT (v • grad)s = — — [r 2 phv r ] -v r -^, (35) 



is calculated to be 



^•« = i$ sW0"-^ (36) 

since the specific enthalpy is generally (i.e., irrespective of the polytropic law) given by 
h = (5/2) (p/p) for an ideal gas. The result is exactly the same as the (effective) Joule 
heating rate, thus confirming a fully advective accretion flow. 

In addition to the heat generation discussed above, the pressure gradient along the flow 
also increases the enthalpy of the fluid. For a unit volume, this amounts to 

"(r, - v r f r = ^M g n K (r ) (^)" 4 sech 4 e = V>, £). (37) 

However, it is more convenient to discuss the energy budget in terms of the quantities 
per unit mass of a fixed fluid element. Then, the Bernoulli sum becomes of our interest. 
Neglecting higher order terms in 3?" 1 in the definition of the sum, we can show that 

K ^_l_GM +h = Q 

2 r 

from our solution. Generally the left-hand side may be a function of position, but in the 
present case it is constant throughout the flow. Therefore, in particular, K is constant 
along the stream lines. This fact also implies that the flow is fully advective. The magnetic 
field plays only a catalytic role in the energy budget. 
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5. INNER AND OUTER EDGES 

The ADAF solution described in the previous section naturally has a finite range of 
validity in the radial direction. The limiting radii for this validity are called the edges. The 
inner edge exists at a radius where the infall velocity becomes comparable with the rotation 
velocity. This means that, at around the inner edge, 9ft becomes of order unity and the 
assumption 5) in §2 becomes invalid. Therefore, we fix the inner edge at the radius where 
3^( r m) = 1- As seen from equation (|52"D, this guarantees the relation a" 1 oc A 2 and that 
3ft(r) > 1 for r > r in . Thus the definition of the inner edge is conceptually clear, but the 
calculation of its explicit expression is generally not so simple. 

On the other hand, the definition of the outer edge is rather vague even in conceptual 
level. There is no general way of defining it. Both the definition of the outer edge and the 
explicit calculation of the inner edge depend on the configurations of external magnetic 
fields. Therefore, we discuss only two typical cases separately. Although similar problems 
have been discussed in K 87, we need to repeat and refine the discussions since the 
r-dependences of the present solution is slightly different from the previous one. 

The first example to be considered is a dipolar external field, because such a 
configuration may well be conceivable around a neutron star or a white dwarf in a mass 
exchanging binary system. In this case, all the advected energy is released finally at the 
surface of the central object. Anyway, from the r~ 3 -dependence of the external dipolar field, 
its strength exceeds that of the induced field in the region within the inner edge. Therefore 
the place at which the boundary value for bg should be fixed is the inner edge, i.e. r = r in , 
and the boundary value is 

B = A (39) 

in 

where \x is the magnitude of a dipole moment and its direction has been assumed to be 
vertical to the midplane of the accretion disk. Substituting this into equation (0), we 



obtain the explicit expression for r in as 
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.2 ,,4 \l/7 



'*^-\ . (40) 

3 GMM 2 V ' 



The position of the outer edge is determined from a requirement of global magnetic-flux 
conservation. Namely, we require that the unperturbed magnetic flux penetrating the 
equatorial plane on the outside of the inner edge is compressed by the presence of an 
accretion flow within the region between the inner and outer edges, i.e., 



B / -f 2vrr dr = / b e 2vrr dr. (41) 

Jr in r Jr in 

In this case, therefore, the outer edge corresponds to a screening radius. The ratio of r out 
to r in is obtained from this equation as 

C = !^ = (1 + A- 1 ) 2 ~A- 2 . (42) 

'in 

The second example, which is of our present interest, is a uniform external field 
perpendicular to the equatorial plane. This is an idealization of the configurations expected 
in various AGNs. In this case, the external field dominates over the induced field at large 
radii. Therefore, the place at which the boundary value for bg should be fixed is the outer 
edge, i.e. r = r out . From equation (|33|) and the relation 9£(r out ) = C^ 2 > we have an implicit 
expression for r out 



/ 3 GMM 2 \ 1/5 , s 

Also in this case, the distance ratio of the outer to inner edges is obtained from a 
consideration of global magnetic-flux conservation. We assume that all the poloidal flux 
generated by the disk and penetrating the disk region should close in the region outside the 
outer edge. Then, the decrease in the flux of the external field in the outer region (r > r out ) 
is balanced by the increase in the inner region (r < r- m ), i.e., 

/•rout „ 

vr in 2 b r {r in ) = / b e 2vrr dr, (44) 
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where the typical field strength in the inner region has been approximated by b r . From this 
equation, we obtain the same result as (§^) again. Therefore, when A 2 <^ 1, equation fl43|) 
reduces to 

(45) 



out 



I B 4 



Here, we briefly return to the energy budget in the disk and discuss it, this time, from 
a global point of view. It has been shown in the previous section that the rate of increase of 
the enthalpy in a unit volume through the current dissipation and the plasma compression 
is 5gj + (r, £). Integrating q^ + over the whole volume of the disk region, we obtain 



Q,i 



dr / d£ 2irAr 2 q+(r, £) 



GMM 
6 



1 



r out 



(46) 



where B has been eliminated by equation (J33| ). The result depends neither on the type 
of the external field discussed above nor on the explicit expressions for r in and r out . As 
expected for a fully advective flow, the increase of the enthalpy in the whole disk region is 
exactly balanced by the net enthalpy flux coming out of the disk region: 



^enth = / hpw ■ dS 



5GMM 
6 



1 



r out 



5Qj 



(47) 



where dS is the surface element of the disk region. 



On the other hand, the kinetic energy flux coming out of the disk region is 



rot 



M 

T 



v v 2 (r- m ) ~v v 2 (r ont ) 



GMM 
6 



r out 



(48) 



Thus, it has turned out that the net input of the gravitational energy in the disk region per 



unit time, -GMM (n 



out 



; ), is converted into the kinetic energy of quasi-Keplerian 



rotation (by 1/6) and the enthalpy which has been added as heat (by 1/6) and also as 
compression (by 2/3). 
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6. SCALING FOR ACTIVE GALACTIC NUCLEI 

Since we are now interested in the AGN activities, the quantities in our model are 
suitably scaled here to such circumstances. The model is completely specified by three 
boundary values and one parameter, i.e., mass of the central black hole M, mass accretion 
rate M, strength of the global magnetic field Bq and the dissipation parameter A. They are 
normalized, respectively, by 1O 8 M , the Eddington accretion rate Me, 1 Gauss and 0.1: 

m "To?V ™ = W E > 6o "§' ^bT (49) 

The explicit expression for the Eddington accretion rate is 

ME=L4xl ° 25 (i^k) gs " (50 > 

Since the reference radius should be fixed at the outer edge, i.e. r = r out , we have 
3?(r out )=A- 1 . 

In terms of the above non-dimensional quantities, the boundary values of other 
quantities appearing in the model are expressed as 

r out = 0.96 x 10 17 bo 4/5 m 2/5 m 3/5 cm, (51) 

Mrout) = 10 r 1 b G, (52) 

M^uO = 3.7 x 10 8 b 2 /5 m- l/5 m 1/5 cm s" 1 , (53) 

p(r-out) = 4.0 5' 2 bl dyne cm" 2 , (54) 

T(r ont ) = 2.7 x 10 8 bf J m~ 2/5 m 2/5 K, (55) 

p(r out ) = 8.8 x 10" 17 cT 2 b 6 /5 m 2/5 m~ 2/5 g cm" 3 . (56) 

All other induced quantities can be calculated from the set of analytic solutions given 
in §4 and their boundary values given above. For example, integrating the density over the 
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disk height, we obtain the surface density which becomes independent of r in our model: 

/oo 
p(r, £)rA d£ = const. = E(r out ), 
-oo 

£(r out ) = 1.7 x 10 cT 1 b 2 J b m 4/5 m 1/5 g cm" 2 . (57) 

The dependence of £ oc M 4 / 5 suggests that the disk is secularly stable (e.g., Kato, Fukue & 
Mineshige 1998). Further integrating S(r) over the radius, we have the total mass in the 
disk as 

M D ~ 7rr 2 t S = 4.9 x 10 35 cT 1 b~ 6/5 m 8/5 m 7/5 g. (58) 

In order for our analytic solution to serve well as a model of ADAF in a global 
magnetic field, the radiation energy which is expected to come out of the disk should be 
really negligible in the physical situations of our interest. This point will be checked here. 
The opacity due to the electron scattering is k cs = 0.4 cm 2 g _1 and the free-free opacity is 
written as kr = 0.6 x 10 23 p T~ 7 ' 2 cm 2 g -1 . In view of the scaled values for p and T, one 
can safely conclude that the former opacity is dominant over the latter. Then, the optical 
depth of the disk becomes r-independent: 

/oo 
«(r, p( r , rA d £ - tcsS(r 2 ) 
-00 

= 3.4 5- 1 b 2 /5 m 4/5 m 1/5 . (59) 

Provided that m <S 1, it can be confirmed from this expression that the disk is everywhere 
optically thin unless its opening angle is extremely small. 

As a representative of radiation losses, we evaluate the Bremsstrahlung from the disk. 
This is because it can be evaluated by a simple formula and other mechanism such as the 
synchrotron loss seems to have a similar order of magnitude (Narayan et al. 1995, Manmoto 
et al. 1997, Narayan et al. 1998). The radiation flux from the disk surface is 

r N -3/2 



/oo 
q B (r, £)rAd£ = F B -(r, 
-00 



out, , 

r out 
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F B ~(r out ) = 2.0 x 10 7 5~ 3 b 2 mm ergs cm" 2 s -1 , (60) 

where ^(r, £) is the volume emissivity. On the other hand, height integrating the rate of 
enthalpy increase, we obtain for energy input rate at a radius r as 

/oo / r \ ~ 3 ' 2 

fej(r, + Mr, £)] r A d£ = Q e + h (r out ) — , 

-oo V'out/ 

g e + h (r out ) = 2.9 x 10 7 6q 2/5 ™~ 1/5 m 1/5 ergs cm" 2 s" 1 . (61) 



Comparing this with the expression (|60|), we can see that the radiation loss is actually 
negligible under the circumstances with m< 1 unless the disk is extremely thin. 

Although our model assumes, for simplicity, a common temperature for both the 
electron and ion components, actually it may take different values for different components 
as is the case in the viscous ADAF models. Such a difference in temperature, however, may 
be less important in the resistive ADAFs because the (effective) Joule dissipation will heat 
mainly the electron component (Bisnovati-Kogan & Lovelace 1997) which is easier to cool 
down by radiation than the ion component. In any case, that is beyond the scope of the 
present paper. 

7. DISCUSSION 

In this section, we discuss main issues which may be raised about our model, 
introducing speculations to some extent. The first is about the accuracy of our set of 
analytic solutions. As described in §4, the separation of variables are performed only 
approximately. It is accurate at the midplane of a disk, but is invalid for large £. The errors 
in the angular dependence may exceed 50% for |£| > 0.5, and certain inconsistencies appear 
beyond tanh |£| = l/v3 where (d 3 /d£ 3 ) tanh£ changes its sign. One may therefore identify 
there (i.e., |£| ~ 1) as the surfaces of a disk. Within this range, the set of solutions serves 
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well for most of our purposes because we are usually interested not in the exact angular 
dependences but in their qualitative behavior. 

Angular momentum is expected to be carried away from the disk by magnetic stresses 
along the externally given poloidal magnetic lines of force. In the case of a dipole type 
external field it is transferred to the central star, and in the case of a uniform external 
field, to distant regions. The rate of angular momentum extraction from the upper and 
lower surfaces of a disk can be calculated from our solution. The result is proportional to 
M\\JGMr ou t — \JGMr\r\ , and the proportionality constant depends on the definition of the 
surfaces £ = ±£ - Such an ambiguity is also a reflection of the approximate nature of the 
solution. 

Although a definite answer cannot be obtained in the present stage of approximation, 
the above picture for angular momentum extraction seems to be physically plausible also 
from a global consideration of the accretion process. A schematic picture is shown in Fig. 
1 for a uniform external field. This field is an idealization of such a large-scale field of 
galactic origin whose sources are in distant regions. The rotation of the accreting plasma in 
this poloidal field acts as a dynamo in a DC-circuit and drives a radial current in the disk. 
We expect this current to close its circuit around distant regions and reaching to the polar 
regions (as for the supply of plasma to this regions, see below). The presence of such a 
globally circulating poloidal current system guarantees the generation of toroidal magnetic 
field outside the disk, which is necessary for the transfer of angular momentum by magnetic 
stresses. 

Thus, the angular momentum can be extracted through a nearly vacuum space, but 
we do not intend to deny the possibility of the extraction also by wind type outflows (e.g., 
Blandford & Payne 1982). Rather, an accretion flow with a divergent component in the 
^-direction has been obtained also in our scheme when the simplifying assumption of vg — 
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is relaxed (Kaburaki 1987). 

The j p x h v force acting on the poloidal current j p whose origin has been discussed 
above contributes to the confinement of accreting plasma towards the equatorial plane and, 
possibly, also to the collimation and acceleration of a bipolar jet system (see, Kaburaki 
& Itoh 1987). The former effect has been described in the present paper as the magnetic 
confinement of an accretion disk (the ^-component of the equation of motion). The latter 
two effects are interpreted as the pinch effect on the return current near the polar axis and 
the acceleration by magnetic pressure-gradient force due to b v , respectively. 

Therefore, the development of a toroidal magnetic field in the middle-latitude regions 
seems to naturally explain such an association of disk-like structure and bipolar-jet structure 
as is frequently observed in many astrophysical situations. Interestingly enough, essentially 
the same story as the above for the formation of a disk-like structure and for the collimation 
and acceleration of a jet applies even to the case of disk-like outflows which are expected 
to exist around rapidly rotating, magnetized stars including pulsars (Kaburaki 1989). This 
is because the sense of the j p x h v force does not change even in that case, in spite of the 
change in the sense of b^ (and therefore that of j p ) which is required to realize an outflow 
(i.e., angular momentum should be added to the disk). Thus, the disk-jet association is 
understood naturally in the same scheme even in the case of outflows. 

The source of plasmas in the polar regions is of course an accretion disk. Within the 
inner edge of an accretion disk, there develops a strong poloidal magnetic field swept by the 
accretion flow. This field may form a core of a bipolar-jet system. The accretion flow would 
be decelerated within the inner edge by the strong magnetic field perpendicular to it and, 
although main fraction of it may be eventually swallowed by the central black hole, some 
part may be turned its direction to follow the core field and then accelerated along it (see 
Fig. 1). 
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It may seem at first rather curious that, as far as our solution is concerned, no energy 
is extracted from the disk associated with the angular momentum extraction by magnetic 
stresses. As confirmed by substituting our solution into the leading-order expressions in A 
of the electric field, 

c 1 (96 1 1 

E r = — KT + -v v be, E e = — (vpbr-Vrbp), E v = 0, (62) 

Aira A r at, c c 

we have E = in the disk and therefore the vanishing Poynting flux. The situation, however, 
is quite analogous to that in the viscous ADAF models in which angular momentum is 
transported radially outwards by viscous stresses but energy is dissipated locally in the flow. 

The vanishing of the Poynting flux is a consequence of idealization that the external 
electric loads in which some part of energy may be consumed is negligible in determining 
the structure of an accretion flow. In such a situation, there is no need to transport energy 
to the outside. This is an opposite extreme of ideal MHD disks. Since no energy dissipation 
can occur in an ideal MHD disk, the whole power input is carried away in terms of the 
Poynting flux to distant regions where the dissipation is assumed to occur. In view of the 
above discussion of global configurations, at least, the jet should be included as an external 
load in a more realistic determination of the current in the disk. However, of course this is 
beyond the scope of the present work. 

The global MHD stability of accretion disks in the resistive ADAF model may be an 
important issue. This is because the ratio of the toroidal to the poloidal magnetic fields 
which is represented by 3? is assumed to be large (3? > 1 and 3? 2 ^> 1) in the model, while 
the configurations with too large 3ft are generally believed to be unstable. Since 3? increases 
from ~ 1 at the inner edge like r 1 / 2 towards the outer edge (3ft(r out ) = A^ 1 in the case of 
uniform external field), the outermost regions are apt to be unstable to the instabilities 
of helical type. Although the determination of the critical value for 3? in actual resistive 
ADAFs is again beyond the scope of this paper, it can be said qualitatively that only not 



-23- 



so widely extended disks and not so thin disks are expected to be maintain safely. 

As for the local MHD stabilities, we have already performed an analysis based on the 
background solution obtained in the present paper and have found that there are actually 
growing modes. Although further investigations including non-linear evolutions are needed 
for a definite conclusion, the instabilities are expected to be not so strong as to destroy the 
whole structure. Rather, they may cause a turbulence in the flow which is the origin of 
some types of anomalous resistivity or viscosity (e.g., Balbus & Hawley 1991). The results 
of the local analysis will be reported elsewhere. 

A. APPENDIX 

Although a reduced Keplerian velocity has appeared in our resistive ADAF model, the 
rotation law may be different for different geometries of the accretion flows. Here, this point 
will be demonstrated by considering two typical cases: one is the disk-like flow of constant 
height (i.e., H = const.) and the other is that of constant opening angle (i.e., A = const.). 
For simplicity, we assume that the disk structure is maintained only by the gravity and that 
the infall velocity is negligibly small compared with the rotation velocity. 

When the disk height is constant, the pressure force is almost vertical to the equatorial 
plane and it is convenient to adopt cylindrical coordinates (w, <p, z). The equations of force 
balance in a poloidal plane are, therefore, 

(Al) 

, (A2) 

r z r p oz 

where r = [w 2 + z 2 ) 1 / 2 and the pressure term in the w-component has been neglected. 
Equation (|A2|) implies a gravitational confinement of the accreting plasma. 



GM w 


= V A 


r 2 r 


zu ' 


GM z 


1 dp 



-24- 



Equation (|A1|) yields the Kepler angular velocity 



/GM\ 1/2 
£l(w, z) = ( — r J = fi K (r). (A3) 

However, the plane of rotation is always parallel to the equatorial plane, so that the center 
of rotation does not coincide with the center of gravitational attraction except for the 
rotation in the disk's midplane. Since the angular velocity is Keplerian, the centrifugal 
force is generally smaller than the gravity at an arbitrary r: 

v% GM . n ^ GM 

-£■ = sin 9 < . (A4) 

w r r 

From equation ( |A2| ) we have a rough estimate for the ratio of sound velocity to the Kepler 
velocity: 

^ ~ ^ • ( A5 ) 

Since if is a constant, this suggests a temperature variation of the form T(r) oc r -3 . 

When the opening angle of a disk is constant, the pressure force is almost along the 
^-direction and it is convenient to adopt spherical polar coordinates (r, 6, if). The force 
balance equations are 

where the pressure term in r-component has been neglected compared with the gravity. 
Equation ( |A7| ) implies a centrifugal confinement of the accreting plasma. 

Equation ( |A6| ) yields the Kepler rotational velocity 

/GM\ 1/2 
«*(*•,*) = (—) =vk(v). (A8) 

The plane of rotation is again parallel to the equatorial plane, so that the center of rotation 
does not coincides with the center of gravitational attraction except for the rotation in the 
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disk's midplane. Since the rotational velocity is Keplerian, the centrifugal force generally 
exceeds the gravity at an arbitrary r: 

< = GM^_ > GM 

zu r sin# r 



From equation (|A7|) we have a rough estimate for the ratio of sound velocity to the Kepler 
velocity: 

C * >° tA| (A10) 



A 



-1 



Since A is a constant, this suggests a temperature variation of the form T(r) oc r 

From the above discussions, we can see that the geometry of disks has essential effects 
on their properties. The rotation law is different between the disks of constant H and 
constant A. The radial dependences of the temperature is also quite different. The same 
dependence as the virial temperature is obtained only in the disks of constant A, but thin 
disks can be realized there only when the force of confinement other than gravity is present 



(i.e., equation (jA7|) cannot hold when A <C 1 even if the temperature is nearly virial). In 
the case of constant H, on the other hand, thin disks are realized if the temperature is 
much smaller than the virial one. 
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Fig. 1. — Schematic drawing of the global geometries of magnetic field and plasma flow. 
A poloidally circulating current system (j p ) driven by the rotational motion of accreting 
plasma generates a toroidal magnetic field b^ in addition to a nearly uniform external field. 
The presence of this toroidal field outside the disk guarantees the magnetic extraction of 
angular momentum from the disk. This field acts to confine the accreting plasma towards 
the equatorial plane and also has a tendency to collimate and accelerate the plasma in the 
polar regions. If the condition is favorable, the plasma in the polar regions may form a set 
of bipolar jets. 



